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Abstract 

We generalize a "one eigenstate" theorem of Levy, Perdew and Sahni 
(LPS) [l] to the case of densities coming from eigenmixture density op- 
erators. The generalization is of a special interest for the radial density 
' functional theory (RDFT) for nuclei [2], a consequence of the rotational 

, invariance of the nuclear Hamiltonian; when nuclear ground states (GSs) 

■ have a finite spin, the RDFT uses eigenmixture density operators to sim- 

_! ' plify predictions of GS energies into one-dimensional, radial calculations. 

, We also study Schrodinger equations governing spin eigendensity matrices. 

, The theorem of Levy, Perdew and Sahni [Tj may be described as follows: 

i) let Ha be a Hamiltonian for A identical particles, with individual mass m, 

- A A 

l>S : Ha = ^[-Ti^A,- /(2m) + u{n)] + ^ vin,f,), (1) 

ii) consider a GS eigenfunction of Ha, il:{ri,(Ti,r2,(J2, ■■■,rA,crA), where cr^ de- 
notes the spin state of the particle with space coordinates r^, 

iii) use a trace of \'ip){il}\ upon all space coordinates but the last one, and upon 
all spins, to define the density, 

p{r) ^ A ^ J dridr2---drA-i\i'{ri,ai,r2,a2, ■■■,rA-i,crA-i,r,aA)f, (2) 

0"1...<TA 

iv) then there exists a local potential Veffif) so that, 
[-h'Ar/i2m) + veffir)] = {Ea - VpM, (3) 
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where the eigenvalue is the difference of the GS energy Ea of the ^-particle 
system and that, Ea-i, of the {A — l)-particle one. 

Can this theorem be generahzed for densities derived from eigenoperators, 
^ J2''^=i^n\4'n){''Pn\, Corresponding to cases where H has several (A/" > 1) 
degenerate GSs V'n? The degeneracy situation is of a wide interest in nuclear 
physics for doubly odd nuclei, the GSs of which often have a finite spin, and, 
if only because of Kramer's degeneracy, for odd nuclei. In particular, because 
of the rotational invariance of the nuclear Hamiltonian, the density operator of 
interest for the RDFT reads, V = J^m I'0Jm)(V'Jm|/(2J+1), where J and M 
are the usual angular momentum numbers of a degenerate magnetic multiplet 
of GSs ipjM- Actually, more generally, it will easily be seen that the argument 
which follows holds for a degenerate multiplet of excited states as well. 

This paper proves the generalization, by closely following the argument used 
for one eigenstate only [1]. Furthermore, there is no need to assume identi- 
cal particles. No symmetry or antisymmetry assumption for eigenfunctions is 
needed. Let pi and at be the momentum and spin operators for the ith particle, 
at position rj. Single out the A-th particle, with its degrees of freedom labelled 
f and a rather than ta and a a ■ For a theorem of maximal generality, with dis- 
tinct masses, one-body and two-body potentials, our Hamiltonian may become, 
'Ha — Ti-A-i +Va + hA, with 

A-1 A-1 

TCa-i = ^[-~h'^ApJ{2mi) +Ui{fi,pi,ai)]+ ^ Vi-i{fi,pi,ai,fj,pj,aj) , 

i—l i>j — l 

A-1 

Va = ^VAjir,'rj,Pj,aj) , Ha ^ ~h'^Ap/{2mA) + UA{r) ■ (4) 
i=i 

The potentials acting upon the first (^—1) particles may be non local and spin 
dependent, but, for a technical reason which will soon become obvious, those 
potentials acting upon the A-th particle in Va and Ha must be strictly local 
and independent of the A-th spin. For notational simplicity, wc choose units so 
that ?i^/(2myi) — 1 from now on. 

As in the one eigenstate case [1] we select situations where there exists a 
representation in which, simultaneously, the Hermitian Hamiltonian Ha and all 
the eigenfunctions '0(ri, ai, f^i, cta) under consideration are real. This reality 
condition does not seem to be restrictive, in view of time reversal invariance. 

Let Ea be a degenerate eigenvalue of Ha- The degeneracy multiplicity being 
larger than 1, select > 2 of the corresponding eigenfunctions ipn, orthonor- 
malized. Their set may be either complete or incomplete in the eigensubspace. 
The density operators, 

2? = ^ IV'n) (V'nl , ^Wn = l, (5) 
n—1 n—1 

with otherwise arbitrary, positive weights w„, are normalized to unity, TrI? = 1, 
in the A-body space. They are eigenoperators of Ha, namely Ha 'D — Ea T). 
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The partial trace of a X> upon the first {A — 1) coordinates and all A spins, 
= X! dfi...drA-i[ipniri,ai,...,fA-i,aA-i,r,a)]'^, (6) 

n— 1 cJi ...CT^ — iO" 

defines a "density" r, normalized so that J df T{r) — 1. Let now ipnra be, in 
the space of the first {A — 1) particles, an auxiliary wave function defined by. 



4'nra[ri,(7i, . . . , 7^4- 1 , (TA- 1 ) = V«(n,cri, ...,rA-i,(JA-i,r,a)/^/T{r) . (7) 

Note that this auxiliary wave function depends on the choice of the weights Wn ■ 
Now the density operator in the space of the first {A — 1) particles, VL = 
X^no- \4'nr(j) w„ {ipurirl, IS normalized, Tr' Vi = 1, where the symbol Tr' means 
integration upon the first {A — I) coordinates and sum upon the first (A — 1) 
spins. Since this normalization of Vp in the [A — l)-particlc space does not 
depend on r, two trivial consequences read, WpTi'VL — and AfTi'Vi — 0. 
More explicitly, this gives, 

T,naai...aA-i I ^^^■■■'^^^-^ ^" (J^n r (^1 , ^1 , . . . , ryl_i , (7^-1 ) X 

'^f 4>nrcriri,ai, . . . , f^-l , (TA- 1 ) =0, (8) 

and 

J dfi...drA-l Wn {[Vr4'nrcriri,<Ti, ...,fA-l,CrA-l)]'^ + 

n (7 (7\. ..(TA — 1 

(i)nra{'ri,(yi,-.-,rA^l,(JA-l) Ar<i>nra{'ri,(7l,-.-,rA~l,(JA-l)} = 0.(9) 

Then one can rewrite the eigenstatc property, {T-La — Ea) V'n = 0, into, 

[Ha-I +VA+hA-EA) <t)nr.^^. (10) 

This also reads, 

\/t {Ha-I + Va + Ua - Ea) (j)nra " {^r \/t) 4>nrcr = 

2 (Vf V7) • (Vf=-0„p^) + V7 {l^f^(i>nra) ■ (H) 

The right-hand side (rhs) of Eq. (fTT|) occurs because the Laplacian, A^?, present 
in Ha, acts also upon the parameter, f, of (f>nra-- This is where the local, multi- 
plicative nature of ua and vaj in the last particle space is used and avoids the 
occurrence of further terms, that would induce a somewhat unwieldy theory. 

Define, for any integrand ^'nfcr, the following expectation value in the first 
{A — l)-particle space, 

((*nro-)) = ^ J dfi...drA-i'^nrcriri,cri, ■■■,rA~i,aA-i) ■ (12) 
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Multiply Eq. PT|) by (pnra and integrate out the first {A — 1) coordinates and 
spins, to obtain, 

{{€f.)) [KTir] + Ea-1 + C/„.(r) + ua{7^ - i?A - M ^W) 
where £^^^"^(7^) is defined fi'om, 

4>nra{ri,Cri, ...,rA-l,CFA-l) 0n r (''l , , . . . , i , flA-l ) , (14) 

and Una-{r) results from, 

4>nra{ri,Oi, . . . , r A^l , <7 A-l)vAj (r, rj , Pj , a j)4>n r ain , CTl , . . . , r A-1 , (T A~l) -il^) 

The square norm of <t>nfcr in the {A — l)-particle space results from Eq. (|12p 
with = 0^. In Eq. 114)) the expectation value of TCa-i is explicited as the 
sum of the GS energy Ea-i of Ha-i and a positive, excitation energy E^%'^{'r). 
From Eq. (fT5|) . the Hartree nature of the potential Unair) is transparent. 

Keeping in mind that, Vr, the density operator 2?' is normalized to unity, 
namely, that X^ncr ^"(('^nf cr)) = 1' multiply Eq. by w„ and perform the 
sum upon n and cr. This gives, 

[U-'-^r) + Ea-1 + U"'-'{r) + ^^(f) - Ea ^ Ap] ^ - 

Sno-"^" [2(((/)„r<j (Vf=-(/)„f?o.))) • (V^Vt^) + (((/)„ ^=-0. (Ap0„ffo-))) ,(16) 

where the "mixture excitation potential" , 

C/-^(f) = ^ u;„ ^^^.^(r-*) , (17) 

na 

is local and positive and the "mixture Hartree-like potential", 

^Hrt^^ ^ ^ ^^^2 ^^^^^ ^ (18) 

no" 

is also local. Because of the frequent dominance of attractive terms in VAj-, 
it may show more negative than positive signs. Then notice that, because of 
Eqs. dS]), and Eq. ^ with = (/iVc/i, the sum in the rhs of Eq. (|16p . 
Ylina {{4>nra C^ffj^nra))), vauishcs. Note also, from Eqs. (HI), and Eq. p2|) 
with = (pAcj), that, again for the rhs of Eq. p6)) . the following equality holds, 

- ^ W„ {{(j)nra {A.p(j)npa) )) X! "^^ {^r't'nra) ' ( Vf (?!)„ f^<j ) )) , (19) 
na ncr 

where, again, the symbol (( )) denotes the trace Tr', an integration upon the 
first {A — 1) coordinates together with summation upon their spins. The rhs of 
this equation, Eq. (dH), defines a positive, local potential, 

U'^'^'if) =J2wn{{ (yrKfa) ' (V,-0„,-,) )) . (20) 
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Finally, according to Eq. the sum of local potentials, W^-^ = JJ'^^'^ + 

jjHrt _|_ jjkm _|_ drivcs a Schrodinger equation for -y/r, 

+ U^ff{r)] v/7(ry = {Ea - Ea-i) VW) ■ (21) 

This is the expected generalization of the LPS theorem. Note that, if the [A— 1) 
particles are not identical, then Ea-i^ the GS energy of 7iyi_i, means here 
the mathematical, absolute lower bound of the operator in all subspaces of 
arbitrary symmetry or lack of symmetry. In practice, however, most cases imply 
symmetries in the (A— l)-space, and Ea-i means the ground state energy under 
such symmetries. 

For nuclear physics, this generalization can be used in two ways: 

i) The first one consists in considering hypernuclei or mesonic nuclei, where the 
A-th particle is indeed distinct. Theoretical calculations with local interactions 
for the distinct particle may be attempted while non local and/or spin dependent 
interactions for the A—1 nucleons are useful, if not mandatory. Then, obviously, 
the density t refers to the hyperon or the meson and, given the neutron and 
proton respective numbers N and Z, wave functions ipn and 0„ „ belong to 
both iV- and Z-antisymmetric subspaces. The energy Ea-\ is the GS energy of 
nucleus {iV, Z}, a fermionic GS energy, rather than the absolute lower bound 
of the mathematical operator Ti^-i in all subspaces. 

ii) The second one consists in setting all A particles to be nucleons, at the cost of 
restricting theoretical models to local interactions. Such models are not without 
interest indeed, although interactions which are spin dependent are certainly 
more realistic. The antisymmetric properties of the functions ■(/'„ are requested 
in both N- and Z-spaces. If the singled out, A-th coordinate is a neutron one, 
the density r defined by Eq. ^ is the usual neutron density, divided by N\ the 
functions <i>nr(y are antisymmetric in the (TV— l)-neutron space and the Z-proton 
space; the energy Ea-\ now means the fermionic GS of nucleus \N — 1, Z}, not 
that absolute, mathematical lower bound of operator TiN-i.z- Conversely, if the 
A-ih. coordinate is a proton one, then, mutatis mutandis, t is the usual proton 
density, divided by Z, and Ea-i is the GS energy of nucleus {N, Z — 1}. 

In both cases, the Hamiltonians to be used are scalars under the rotation 
group, and, therefore P], the density operators V considered by the RDFT are 
also scalars. Hence, all calculations defining V^^^ reduce to calculations with a 
radial variable r only. 

We shall now extend our previous results to the case where we allow spin 
dependence for all interactions, a most useful feature if all A particles are nucle- 
ons. Polarized eigenmixtures are also interesting and need also be considered. 
Hence, a generalisation of our approach, which uses the "spin-density matrix" 
(SDM) formalism [3] [4], is in order. The Hamiltonian may become, 

A A 

'Yy-'h^^TiK'^mi) + Ui{ri,Si)]+ ^ Vtj{n,ai,fj,aj) . (22) 

1=1 i>j = i 

It allows subtle differences between neutrons and protons, besides the Coulomb 
interactions between protons. More explicitly, there can be two distinct one- 
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body potentials, Un,Up, namely one for neutrons and one for protons, but within 
the neutron space the function Un{ri,ai) obviously will not read Uni{fi, (Ji). Simi- 
larly in the proton space, the Hamiltonian contains terms Up{fi, (T,) rather than 
The same subtlety allows terms Vpp{n,ai,rj,cTj), Vpn{ri,(Ji,rj,aj), 
) and Vpp{'ri,ai,fj,aj), rather than Vppij(fi,(Ti,rj,aj), ... etc. 



(Of course, Vp 



.) But non localities of potentials and interactions, in the 



sense of explicit dependences upon momenta pi, remain absent. 

Then the A-th. particle is again singled out, with degrees of freedom again 
labelled r and a, and the Hamiltonian is split as a sum, K-a-i + + kA, some- 
what similar to the split described by Eqs. For simplicity, we shall use short 
notations, /C, W and fc, rather than /Ca-i, Wa and kA- With two spin states, 
a — ±1/2, for the A-th nucleon, we represent eigenstates ipn as column vec- 

, and operators as matrices such as, W = yy^^ VV^^ 

f ^ and u = ^'^++ . The matrix, K, — ^ ^ 

K__J \_UA h "^A J [0 /C 

scalar in spin space, since K. does not act upon the A-th particle. 

The spin density matrix, p„, results from an integration and spin sum over 

the [A — l)-space of the matrix, -(/;„ 
transposition. 



tors, 

'k. 



tpn , where the superscript denotes 



Pnir) = (( 



X [V'„+ V'n-] )> = (( 



(V^n+) 



)) • (23) 



It depends on the last coordinate, r, and its matrix elements are labelled by 
two values, {a, cr'}, of the last spin. For an eigenmixture one defines, obviously, 
^(r) = Sn^nPn(^); and the trace in the last spin space, [0++(r) + 0__(r)], is 
that density, r(r), defined by Eq. ([6]). 

The SDM, 0, is symmetric and positive semidefinite, Vr. Except for marginal 
situations, it is also invertible, in which case there exists a unique inverse square 
root, also symmetric and positive. Define, therefore, a column vector </)„^ of 
states in the {A — l)-space according to, 



(24) 



Then the following property. 



5]«;„((0_^x^_^))=^-^(r) ^«;„((V„xV'„))U-^(r) = i, (25) 



holds Vr. Here 1 denotes the identity matrix. Hence, the following gradient and 
Laplacian properties also hold, 



^ w„ (( ( Vf 

n 



t^7i r — 



'n r + 

'nr — 



X [(t>nr+ 4>nr- ] )) + 
x(VF[0„r+ 0«r-]))) = 0, Vf, (26) 
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and 



f — 



f — 



{yp[(t)„p+ (t)„p^ ]) )) + 

x{Ap[(t>nr+ 0, Vf . 

The eigenvector property, (/C + W + fc — i?^ l) = 0, also reads, 



(27) 



^5 (f) . (V^. (/.„,^) + 9i (f) (A,- . 



(28) 



Right-multiply Eq. (I28p by the row vector, integrate and sum over the 

{A — l)-space, weigh the result by Wn and sum upon n. Because of Eq. (f25|) . 
the weighted sum of averages over the {A — l)-space simplifies into. 



U (f) + {Ea-1 -Ea)1+U (f) + u{r) - Ap\ 0^ (f) = E„ x 
{2[v.^0"^(f)] • {{{\/p^^p) x^l^)) +eHr) {{{Ap^^p) xt^,-))} , (29) 



with 



and 



l^'''"'ir) = Wn (( (/C 0„^=.) X (/)„^. )) - Ea-1 1 , 



(30) 



(31) 



The rhs of Eq. (|29p can be simplified, but less than that of Eq. ([T5)) . Indeed 
Eq. (pS)) does not imply that the coeflicient of Vf ^2(r) vanishes. In fact this 
coefficient is. 



(((Vf0„p+)0„f_)) 
(((VF0„f?_) (?!)„f-+)) 



(32) 

and Eq. (|26p shows that the matrix on the rhs is antisymmetric. In turn, from 
Eq. ([27|l . the "Laplacian induced coefficient" in the rhs of Eq. (|29p may be 
hsted as, 



T,n^n{{{Af<l>nr) X 0„f.)) 



1^71 r + 



i>nr- 



X [ (j)nr+ 4>nf- 



X {Ap[4>n?+ (jjnr- ]) 



»(33) 
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In the rhs of Eq. the similarity of its first term with potential Jj'^™, Eq. 

Ijl9p . is transparent. Also transparent is the antisymmetry of the second term. 
With the definitions, 



and 



(( (Vf='0„f?+)0„f?_ )) 



(34) 

(35) 

(36) 



the Schr5dinger equation for the square root of the spin density matrix reads, 



—exc —Hrt _ —kin —ajit . 

U +U +U+U -~U -A 



>i~U'"'^-V9i ^ {Ea~Ea-i)0^. (37) 



This concludes our generalizations of the LPS theorem. On the one hand, 
see Eq. (HU, we obtained for the square root density of an eigenmixture a 
Schrodinger equation, most similar to the LPS equation. On the other hand, 
at the cost of a slightly less simple result, we also obtained, see Eq. ([37|) . an 
LPS-like equation that drives the square root of the spin density matrix. 

This opens a completely new zoology of local, effective potentials, of which 
very little is known, but the interest of which is obvious, since they drive a 
reasonably easily measurable observable, the square root of an eigenmixture 
density, which depends on one degree of freedom r only. A connection of such 
potentials with optical potentials, or rather their real parts, is likely, but yet 
remains an open problem. 
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